Theory for atomic diffusion on fixed and deformable 
crystal lattices 



Eliot Fried 

Department of Theoretical and Applied Mechanics 
University of Illinois at Urbana- Champaign 
Urbana, IL 61801-2935, USA 

Shaun Sellers 

School of Mathematics 
University of East Anglia 
Norwich NR4 IT J, England 

Abstract. We develop a theoretical framework for the diffusion of a single un- 
constrained species of atoms on a crystal lattice that provides a generalization of 
the classical theories of atomic diffusion and diffusion-induced phase separation 
to account for constitutive nonlinearities, external forces, and the deformation of 
the lattice. In this framework, we regard atomic diffusion as a microscopic process 
described by two independent kinematic variables: (?) the atomic flux, which reckons 
the local motion of atoms relative to the motion of the underlying lattice, and (ii) the 
time-rate of the atomic density, which encompasses nonlocal interactions between 
migrating atoms and characterizes the kinematics of phase separation. We introduce 
generalized forces power-conjugate to each of these rates and require that these forces 
satisfy ancillary microbalances distinct from the conventional balance involving the 
forces that expend power over the rate at which the lattice deforms. A mechanical 
version of the second law, which takes the form of an energy imbalance accounting 
for all power expenditures (including those due to the atomic diffusion and phase 
separation), is used to derive restrictions on the constitutive equations. With these 
restrictions, the microbalance involving the forces conjugate to the atomic flux pro- 
vides a generalization of the usual constitutive relation between the atomic flux and 
the gradient of the diffusion potential, a relation that in conjunction with the atomic 
balance yields a generalized Cahn-Hilliard equation. 

Dedicated to Roger Fosdick, whose fundamental and broad-reaching 
work in continuum thermomechanics has provided us with abundant 
inspiration. 



1. Introduction 



The diffusion of a single unconstrained atomic species on a fixed crystal 
lattice is commonly modeled by the Cahn-Hilliard equation (see, for 
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example, Cahn |T], §, | and Cahn & Hilliard |, |, §) 

z> = Div Grad i/ - D 2 (u) Grad (Div Grad u)) , (1.1) 

with ^ the atomic density, D\ the nonnegative density-dependent dz/- 
fusivity, and -D2 a small, nonnegative density-dependent coefficient. 
The first term on the right side is due to simple diffusion, whereas 
the second one is a higher order term that takes into account spatial 
inhomogeneities on a fine scale, such as those that may arise during 
processes involving phase separation. This equation can be derived by 
adjoining to the atomic balance 

v = -Div j, (1.2) 

the constitutive relation 



j = —k(u) Grad fi, (1-3) 

determining the atomic flux as proportional to the gradient Grad \x of 
the diffusion potential fi through a nonnegative scalar mobility k that, 
in general, depends on the atomic density. In simple atomic diffusion, 
the diffusion potential is given by the derivative of a constitutive re- 
sponse function ip that determines the energy density ip in terms of the 
atomic density, viz., 

H = (1.4) 
in which case the diffusive atomic flux takes the form 



j = —D{v) Grad v, (1.5) 

with diffusivity D determined by 

D(v) = n{u){jj"{v). (1.6) 

However, when strong spatial inhomogeneities in the density may be 
present, the response function determining the energy density is taken 
to depend on both the density v and its gradient Grad v, so that the 
diffusion potential is given by the variational derivative 

dip dip 
/,= -(,, Grad v) - Div ( ^dv) ^ Gmd V » ' (L7) 

Often, in this context, it is presumed that ip has the particular form 
ip(u, Grad v) = f(y) + \a \Grad v\ 2 , (1.8) 
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with a a positive constant that is usually interpreted as proportional 
to a squared characteristic length of interaction and, hence, taken to 
be small, so that the atomic flux becomes 

j = -L>iO) Grad v + D 2 {y) Grad (Div Grad u), (1.9) 

with coefficients 

Dx(y) = K(v)f"(y) and D 2 {y) = an{v). (1.10) 

Recently, Gurtin M presented a thermodynamically consistent frame- 



work that provides a broad generalization of (1.1) that accounts for 
constitutive nonlinear ities, kinetics, and deformation of the underlying 
lattice. This framework is based on the belief that to each independent 
kinematic process there should correspond a system of power-conjugate 
forces subject to a generalized force balance.^] Specifically, Gurtin |?j 
identifies atomic diffusion and lattice deformation as distinct kinemat- 
ical processes, with the density-rate v being a generalized velocity 
associated with atomic diffusion and the deformation-rate y being the 
velocity of the underlying lattice. Associated with v is a system of 
generalized forces that satisfy a scalar microforce balance. This balance 
is a postulate distinct from and ancillary to the conventional force 
balance involving forces that act power conjugate to the velocity of 
the underlying lattice. With appropriate thermodynamically consistent 
constitutive equations, the microforce balance yields an expression for 
the diffusion potential that generalizes (p.. 7]). Finally, the generalized 
Cahn-Hilliard equation arises on combining with the atomic balance 
the thermodynamically derived expression determining the atomic flux 
as proportional to the gradient of the diffusion potential. 

Although this approach provides a powerful and appealing means 



of deriving (LI) and various generalizations thereof, we believe that it 
is lacking in the sense that it accounts only for one aspect of the kine- 
matics associated with atomic diffusion. In particular, we take the view 
that the diffusive flux of atoms, as described by j, should be considered 
as a kinematic process distinct from that described by the density-rate 
v. Specifically, we conceive of atomic diffusion as a microscopic process 
described by 

(i) the atomic flux j, which measures the local motion of atoms rel- 
ative to macroscopic motion of the lattice as characterized by the 
velocity y, and by 



1 A foundation for this viewpoint can be devised using the principle of virtual 
power, see, for example, Germain Antman & Osborn Maugin jic|] , and 
Fremond Jill. We follow here an alternative, but essentially equivalent, approach 
taken by Ericksen ]l2|] ^Goodman & Cowin |H| , Capriz & Podio-Guidugli |l^] , Capriz 
||, Fried & Gurtin (fij, 0, Fried |@, and Gurtin |]. 
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(ii) the rate v of the atomic density, which encompasses nonlocal inter- 
actions between diffusing atoms and characterizes the kinematics 
of phase separation. 



Adopting the view promoted by Fried & Gurtin 16, 17], Fried jig], and 
Gurtin |?J, we arrive at an alternative derivation of the Cahn-Hilliard 
equation (|P]), where, to account for power expenditures associated 
with the microscopic motion of atoms, we introduce two separate sys- 
tems of forces — one being power-conjugate to the atomic flux j and 
the other being power-conjugate to the rate v of the atomic density. 
Further, we require that these forces be consistent with separate micro- 
force balances, one vectorial and involving the forces power-conjugate 
to j and the other scalar and involving the forces power-conjugate to 

We first focus on the development of a simple theory, where, to 
highlight the roles of the vectorial and scalar microforces, the lattice 
is fixed. Next, we extend the theory to include deformation of the 
lattice. For simplicity, we restrict our attention to processes involving 
only atomic diffusion and lattice deformation, ignoring thermal and 
other effects, so that the first and second laws of thermodynamics are 
replaced by an energy imbalance that accounts for all relevant power 
expenditures. As a further simplification, we ignore all forms of inertia. 
We employ a purely referential description and consistently use direct 
notation (see, for example, Gurtin J2(J). 



2. Simple theory for atomic diffusion on fixed lattice 

2.1. Balance and imbalance laws 

The theory is based upon the following laws: 

• atomic balance; 

• flux-conjugate microforce balance; 

• density-rate-conjugate microforce balance; 

• energy imbalance. 

In contrast to standard formulations of diffusion in rigid bodies, we 
include in this list of balances two microforce balances, which will be 

2 Previously, Fried & Sellers developed a theory for solute transport in which 
the diffusive mass flux is interpreted as a generalized velocity and corresponding 
generalized forces subject to an ancillary force balance are introduced. 
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seen to provide appropriate generalizations of commonly assumed con- 
stitutive relations. In particular, the flux-conjugate microforce balance 
provides the generalization of ( |1.3| ), whereas the density-rate-conjugate 
microforce balance provides the generalization of (1.7). 



We articulate our basic laws of balance and imbalance in global 
form over a generic subregion V of the region B occupied by the body 
in a fixed reference configuration. The outward unit-normal field to dV, 
directed outward from V, is designated by n. 



2.1.1. Atomic balance 
We introduce the fields 



v atomic density, 

j atomic flux, 

m external atomic supply, 



in which case the integrals 



v dv, J j-nda, and J m dv (2-1) 

V dV V 

represent, respectively, the number of atoms in V , the number of atoms, 
per unit time, added to V by diffusion across dV, and the number of 
atoms added, per unit time, to V by agencies external to 

Atomic balance is the postulate that, for each region V and each 
instant, the rate at which the number of atoms within V changes with 
respect to time be equal to the rate at which atoms are added to V by 
the diffusion across dV and by external supplies to V: 



J v dv= — J j n da + J m dv. (2-2) 

V &P V 

Using the transport and divergence theorems, we obtain the equivalent 
local field equation 

i) = — Divj + m (2-3) 

enforcing atomic balance. 

3 The supply m should not be confused with a reaction term. Rather, we view 
m as an external supply that accounts for the insertion of material by an external 
observer. Its role is analogous to that of the external heat supply that is routinely 
included in modern continuum-thermodynamical statements of energy balance (see, 
for example, Truesdell ^J). By introducing the external mass supply m, we follow 
Gurtin 0. 
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2.1.2. Flux- conjugate microforce balance 

Associated with the evolution of the atomic flux j, we introduce a 
system of power-conjugate microforces, consisting of 

£ flux- conjugate microstress, 

h flux- conjugate internal body microforce, 

I flux- conjugate external body microforce, 

in which case the integrals 

J En da, j h dv, and J £dv (2.4) 
TV V 

represent the flux-conjugate microforces exerted on the region V by the 
flux-conjugate microtraction distributed over dV, by agencies within V, 
and by flux-conjugate agencies external to V , respectively. 

Flux-conjugate microforce balance is the postulate that, for each 
region V and each instant, the resultant, taking both internal and 
external sources into account, of the flux-conjugate microforces acting 
on V vanishes: 

J Snda + J (h + t) dv = 0. (2.5) 

ar v 

Using the divergence theorem, we arrive at the equivalent local field 
equation 

Div S + h + £ = (2.6) 
enforcing the flux-conjugate microforce balance. 

2.1.3. Density-rate-conjugate microforce balance 

Associated with the evolution of the atomic density v, we introduce a 
system of power-conjugate microforces, consisting of 

£ density-rate-conjugate microstress vector, 

7r density-rate-conjugate internal body microforce, 

7 density-rate-conjugate external body microforce, 

in which case the integrals 

J ^-n da, j it dv, and J \ dv (2.7) 
V V V 
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represent the density-rate-conjugate microforces exerted on the region 
V by the density-rate-conjugate microtraction distributed over dV, by 
agencies within V, and by density-rate-conjugate agencies external to 
V , respectively. 

Density-rate-conjugate microforce balance is the postulate that, for 
each region V and each instant, the resultant, taking both internal and 
external sources into account, of the density-rate-conjugate microforces 
acting on V vanishes: 

J £-nda+ J(tt + j) dv = 0. (2.8) 
dV V 

Using the divergence theorem, we arrive at the equivalent local field 
equation 

Div£ + 7r + 7 = (2.9) 
enforcing the density-rate-conjugate microforce balance. 



2.1.4. Energy imbalance 

To formulate this imbalance, we introduce fields 

ip free energy density, 
(i diffusion potential, 



in which case the integrals 

J ip dv and J fim dv (2.10) 
V V 

represent the free energy in V and rate at which energy is added to V 
through the external supply of atoms to V. 
Further, the integrals 

J (£ T j + &)-nda and J (£-j + 7 z>) dv (2.11) 
&P V 

provide an accounting of the power expended on V by the tractions 
acting on dV and by the agencies acting external to that region. 
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Energy imbalance is the postulate that, for each region V and each 
instant :[] 



Jipdv< J (E T j + £i>)-n da + J (i-j + 71/ + fj,m) dv. (2.12) 
V dV V 

Using the transport and divergence theorem as before, we arrive at the 
local inequality 

+ tx)v + i- Grad i> - hj + (S + ^1) ■ Grad j > (2.13) 
enforcing energy imbalance. 

2.2. Constitutive equations 

We assume that the free energy density ip, the flux-congugate mi- 
crostress S, the flux-congugate internal body microforce h, the density- 
rate-conjugate microstress and the density-rate-conjugate internal 
scalar body microforce it are determined constitutively in terms of 
the atomic density v, the gradient Grad v of the atomic density, the 
diffusion potential fi, the gradient Grad /j, of the diffusion potential, and 
the atomic flux j: 



(2.14) 



We emphasize that, in contrast to the standard approach to the theory 
of atomic diffusion, neither the diffusion potential [i nor the atomic flux 
j is given by a constitutive relation. 





= ip(v 


Grad v 


H, Grad /x, j), 


£ 


= S^v, Grad v 


Gradfi,j), 


h 


= h{v, 


Grad v. 


fi, Grad > 




= 


Grad v, 


H, Grad fi, j), 


7T 


= 7r(v, 


Grad z/, 


(j,, Grad fx, 3). 



Inserting the constitutive relations (2.14) into the local energy im- 
balance ( 2.13j ) and writing, for brevity, 



{y, Grad v, fi, Grad fi,j), 



(2.15) 



4 To account for the diffusive flux of energy, Gurtin |Q| included in his statement 
of energy imbalance the term 

— J f-ij-n da. 

av 

Since, in our context, the diffusive flux of energy can be interpreted as the surface 
power of the atomic diffusion, our theory accounts for this effect through the term 
with integrand S T j n. 
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we arrive at the functional inequality 

(g W - , + *(,))* + (gjfgfjjM - fw) - o*»* + 

+ ^^^r (2!)-Grad/x+^(z)-j + fa(z)-j-(g(z) + Aii)-Gradj < 0. 

(2.16) 



Hence, following the procedure founded by Coleman & Noll [22] in their 
incorporation of the second law into continuum thermomechanics, we 
find that: 

(i) the constitutive response functions tp and £ delivering the free en- 
ergy density ip and the density-rate-conjugate microstress £ must 
be independent of the diffusion potential fi, the gradient Grad \x 
of the diffusion potential, and the atomic flux j, and obey 

£{y, Grad v) = ^ - (u, Grad v); (2.17) 

(ii) the constitutive response function £ delivering the flux-conjugate 
microstress £ must be independent of the atomic density u, the 
gradient Grad v of the atomic density, the gradient Grad \i of the 
diffusion potential, and the atomic flux j, and obey 

S=-fil; (2.18) 

(Hi) the constitutive response function tt delivering the density-rate- 
conjugate internal body microforce tt must be independent of the 
gradient Grad \x of the diffusion potential and the atomic flux j, 
and obey 

dip 

7r(u, Grad u, /x) = \i — ~q~( v i Grad v); (2.19) 

(iv) the constitutive response function h for the flux-conjugate internal 
body microforce h must be consistent with the residual inequality 

h(u, Grad u, /j, Grad fi,j)-j < 0. (2.20) 

Further, granted smoothness of the response function h, a result due 
to Gurtin & Voorhees p3] yields a general solution 



h = —Z(v, Grad v, fj,, Grad fx, j)j (2.21) 
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of ( |2.20| ), where the reciprocal mobility tensor Z must obey 

3-Z(y, Gradv,/j, Grad > 0. (2.22) 

The governing equations that arise on substituting the foregoing 
thermodynamically consistent constitutive equations in the local field 
equations (^) and ( |2.6| ) expressing flux-conjugate microforce balance 
and density-rate-conjugate microforce balance read 

Z(v, Grad u, fi, Grad fJ,,j)j+ Grad \x = £, 

-(,, Grad v) - Div (^^, Grad v))-» = r J 

Together with the atomic balance ( |2.3| ), ( |2.23| ) form the final govern- 
ing equations of our theory for the diffusion of a single unconstrained 
atomic species on a fixed lattice. 

The foregoing results show that the behavior of a medium of the 
sort considered here is completely determined by the provision of two 
constitutive response functions: 

• ip determining the free energy density as a function of the density 
v and its gradient Grad v; and 

• the reciprocal mobility tensor Z, which, in general, may depend 
on v, Gradv, fi, Grad fi, and j. 

Provided that the external supplies m, and 7 all vanish, that 
the constitutive function ip is quadratic in the gradient of the atomic 
density and twice differentiable with respect to the atomic density, 
and that Z = Ql with £ depending at most on the atomic density, 
the system formed by combining ( 2.23| ) with the atomic balance (2.3) 



yields precisely the Cahn-Hilliard equation (1.1). If, furthermore, ip is 



independent of the gradient of the atomic density, we obtain 



■Gradv, (2.24) 



for the atomic flux j in terms of the density v and its gradient Grad v. 
Comparing ( 2.24| ) with the conventional relation (|1.5|), we obtain an 



expression 



«W = |^ (2.25) 
for the diffusivity. Hence, within our framework, a constitutive relation 



of the form (1.5) can, therefore, be interpreted as an expression of flux- 



conjugate microforce balance, granted that the flux-conjugate internal 
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body microforce is isotropic and linear in the atomic flux and that 
the flux-conjugate external body microforce vanishes. If, moreover, we 
assume that •if) has the particular form 

$(v) = k B 9v\ogv (2.26) 

corresponding to the classical entropic contribution to the free energy 
of a dilute mixture, with /cb Boltzmann's constant and 9 the absolute 
temperature, then (|2.25| ) reduces to 

D(u) = M, (2.27) 

CH 

the standard relation between diffusivity and mobility (see, for exam- 



ple, Nernst [24| and Einstein pql ). More generally, granted that the 
reciprocal mobility tensor Z is invertible, our theory in conjunction 
with the assumption ( 2.26|) yields the diffusivity as a tensor 



D{is, Grad is, fj,, Grad fx, j) = k^OZ 1 (v, Grad is, fx, Grad fi, j). (2.28) 



3. Theory accounting for both atomic diffusion and lattice 

deformation 

3.1. Balance and imbalance laws 



In addition to the mass balance ([2.2]), flux-conjugate microforce balance 
([2.5[), and density-rate-conjugate microforce balance ( |2.8[ ), the theory 
is based upon the following laws: 

• lattice- velocity-conjugate force balance; 

• lattice- velocity-conjugate torque balance; 

• energy imbalance. 

Again, the standard formulation of diffusion in a deformable body is 
supplemented with two microforce balances, which will provide appro- 



priate generalizations of (|1.3j) and (1.7). 



3.1.1. Lattice-velocity-conjugate force balance and 

lattice-velocity-conjugate torque balance 
The motion of the lattice is described by a mapping y of points x in B 
and times t into points y(x, t) of space. We require that Det(Grad y) > 
0. Associated with this motion, we introduce a system of power-conjugate 
forces consisting of 
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S lattice-velocity-conjugate (Piola) stress, 
b lattice-velocity-conjugate body force, 

in which case the integrals 

J Sn da and J b dv (3.1) 

V V 

represent the lattice- velocity-conjugate forces exerted on the region V 
by the lattice-velocity-conjugate traction distributed over dV and by 
lattice- velocity-conjugate agencies external to V , respectively. Further, 
the integrals 

J yxSnda and J yxbdv (3.2) 
V V 

represent the lattice- velocity-conjugate torques exerted on the region 
V by the lattice- velocity-conjugate traction distributed over dV and by 
lattice- velocity-conjugate agencies external to V, respectively. 

Lattice-velocity- conjugate force balance is the postulate that, for each 
region V and each instant, the resultant of the lattice- velocity-conjugate 
forces acting on V vanishes: 

J Snda + J bdv=0. (3.3) 

dV V 

Further, lattice-velocity- conjugate torque balance is the postulate that, 
for each region V and each instant, the resultant of the torques associ- 
ated with the lattice-velocity-conjugate forces acting on V vanishes: 

J yxSnda + J yxbdv =0. (3.4) 

dV V 

In writing fl3.4D , we have assumed that the diffusing atoms do not in- 
duce torques, as would be expected at least for approximately spherical 
atoms. Using the divergence theorem, we arrive at the equivalent local 
field equations 

DivS + b = and (Grad y) S T = S (Grad y) T . (3.5) 
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3.1.2. Energy imbalance 

To account for power expenditures associated with the motion of the 
lattice, we modify the energy imbalance ( p. 12 ) to read 



J ijj dv < J (S T j + £z> + S T y) -nda + J (I ■ j + nv + fj,m + by) dv. 
V dV V 

(3.6) 

Proceeding as before, we arrive at the local inequality 

— ip + (fi — 7r)i> + £-Grad v — h j 

+ S-Grady + (£ + fil)-Gradj > 0. (3.7) 

3.2. Constitutive equations 

We add the lattice- velocity-conjugate stress and the lattice- motion gra- 
dient Grady to the lists of dependent and independent constitutive 
variables, respectively, so that 

ip = ^(v, Grad v, fi, Grad Grad y), 
£ = £{v, Grad v, /x, Grad [i, j, Grad y) , 

h = h(i/, Grad v, fi, Grad fi, j, Grady), ■ 

} (3.8) 

£ = £(z^, Grad v, n, Grad fi,j, Grad y), 
ir = -k(v, Grad v, Grad fi, j, Grad y), 
S = S(v, Grad v, /i, Grad [i, j, Grad y). , 

Inserting the constitutive relations fl3.8|) into the local energy imbal- 
ance ([^7]) and writing 

z = {y, Gradv,fj,, Grad/j,,j, Grady), (3-9) 

we arrive at the functional inequality 

( £f w _ „ + + ( _J*_ w _ j W) . 0rad , + | (2)(i 

Qip dip 
+ W^)} Z) ' Gradfl+[ d{Grady) ){z) ~ S ^)- Grad V 

+ ^(z)-3 + h(z)-j-{±{z)+ l il)-Grad3 < 0. (3.10) 
Hence, proceeding as before, we obtain the results 
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(i) the constitutive response functions ip, S, and £ delivering the 
energy density ip, the lattice- velocity-conjugate stress S, and the 
density-rate-conjugate microstress £ must be independent of the 
diffusion potential /x, the gradient Grad /x of the diffusion potential, 
and the atomic flux j, and obey 



dip 

S(u, Gradv, Grady) = — — — ; — r(u, Gradv, Grady), 

o(Grady) 

dip 

£(v, Grad v, Grad y) = — — — - — -(v, Grad u, Grad y); 

d(Gradu) 



(3.11) 



(ii) the constitutive response function £ delivering the flux-conjugate 
microstress 2 must be independent of the atomic density v, the 
gradient Grad v of the atomic density, the gradient Grad fi of the 
diffusion potential, and the atomic flux j, and obey 

S = -/j,l; (3.12) 

(Hi) the constitutive response function tt delivering the density-rate- 
conjugate internal body microforce tt must be independent of the 
gradient Grad \x of the diffusion potential and the atomic flux j, 
and obey 

dip 

n(u, Grad u,fi) = // — ~q~( v i Grad v); (3.13) 

(iv) the constitutive response function h for the flux-conjugate internal 
body microforce h must be consistent with the residual inequality 

h(u, Grad v, ix, Gradfi,j, Grady)-j < 0. (3-14) 
Analogous to the result ( 2.21| ), we have the general solution 



h = — Z(v, Grad v, fj,, Grad fi,j, Grad y) j (3.15) 

of ( |3.14| ), where the reciprocal mobility tensor Z must obey 

j-Z(u, Grad v, \i, Grad //, j)j > 0. (3.16) 

Further, the constitutive relations are restricted by objectivity; that 
is, under the observer transformation 



v v, Grad v \— > Grad v, j i— > j, 

fx i— > fx, Grad fx i— > Grad fx, Grad y t— > Q(Grad y), 



(3.17) 
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with Q orthogonal, we must have 

S^QS, h^h, £^£ 

TT^TT, 7^7,} (3.18) 

S i— > QS, b ^ b, ip i— > ip.. 

Granted the foregoing thermodynamic results, these transformations 
require, in particular, that 



ijj = ijj(v, Grad v, C), 

Z = Z(u, Grad u, fx, Grad fj,,j, C), 

S = (Grady)^(v, Gradu,C), 



(3.19) 



with C = (Grad y) J Grady the right Cauchy-Green tensor. 

The governing equations that arise on substituting the foregoing 
thermodynamically consistent constitutive equations in the local field 
equations ([T^)i, ( ^9j ) and (|2.6|) expressing lattice- velocity-conjugate 
force balance, flux-conjugate microforce balance, and density-rate-conju- 
gate microforce balance read 

Div ((Grad y)^(v, Grad v, C))+b = 0, 

Z(u, Grad u : //, Grad C) j + Grad [i = 

dip dip 

-(v, Grad u, C) - Div (^^y {», Grad u, C)) - M = 7 . 

(3.20) 

Together with the atomic balance ( |2.3| ), ( |3.20| ) form the final govern- 
ing equations of our theory for the diffusion of a single unconstrained 
atomic species on a deformable lattice lattice. 

We find that the behavior of a medium of the sort considered here 
is still completely determined by the provision of two constitutive re- 
sponse functions: 

• ip determining the energy density as a function of the density u, 
the gradient Grad v of the atomic density, and right Cauchy-Green 
tensor C associated with the lattice deformation; and 

• the reciprocal mobility tensor Z, which, in general, may depend on 
the atomic density v, the gradient Grad v of the atomic density, the 
diffusion potential (jl, the gradient Grad /i of the diffusion potential, 
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the atomic flux j, and the right Cauchy-Green tensor C associated 
with the lattice deformation. 



4. Discussion 

Our theory is based upon introducing generalized velocities for diffusion 
and corresponding power-conjugate forces, which we call microforces, 
into an otherwise conventional continuum-mechanical description. In 
particular, we identify two distinct generalized velocities — a vector- 
valued diffusive atomic flux that characterizes simple diffusion and a 
scalar- valued atomic density rate that characterizes phase separation. A 
systematic, thermodynamically consistent derivation generates a final 
system of governing equations — ( |2.3| ) and ( |2.23| ) for a fixed lattice, and 
( p. 3[ ) and ( |3.20D for a deformable lattice — that constitutes a generaliza- 
tion of common diffusion equations, a generalization that accounts for 
constitutive nonlinearities and external forces acting on the diffusing 
atoms. 

We find that the constitutive response of the material is determined 
by the provision of two functions, one for the energy density and the 
other for the reciprocal mobility tensor. In particular, the Cahn-Hilliard 
equation ( |1 . 1[ ) results on assuming that the energy response function 
is quadratic in Grad v and that the reciprocal mobility tensor is pro- 
portional to the identity tensor by a scalar coefficient that is at most a 
function of v. Importantly, we arrive at our theory without introducing 
a constitutive relation such as (|l.3| ) for the diffusive flux. 

Additionally, our proposed format also leads to a generalization, 
( [2.25 ), of the classical relation between diffusion coefficient and viscos- 



ity. We find that the classical form for the relation ( |2.27| ) holds when 
the free energy is given by the special constitutive relation ( |2,26| ). 

A distinction between the final governing equations of our theory 
and those arising in more conventional derivations is that, in our ap- 
proach, a vector- valued microforce balance provides a generalization 
between diffusive atomic flux and gradient of the diffusion potential. 
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